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Abstract — Under simulated physiological conditions 
(pH=7.40), the interaction between Mobic and lysozyme was 
studied by synchronous fluorescence spectroscopy, 
ultraviolet spectroscopy, circular dichroism spectroscopy 
and molecular docking simulation technique.The results of 
spectroscopy showed that the fluorescence of lysozyme was 
statically quenched by Mobic, the number of binding sites 
was about 1, and the conformation of lysozyme was 
changed.The thermodynamic parameters obtained from the 
van't Hoff equation show that the Gibbs free energy A G<0 
showed that the reaction between them was spontaneous, 
and AH<0, AS>0, indicating that the hydrophobic force 
plays an important role in the formation of Mobi c-lysozyme 
complex.The results of molecular docking showed that the 
binding site of Mobic was close to the active site composed 
of Asp48 and Glu35 residues, indicating that Mobic could 
change the microenvironment of amino acid residues at the 
catalytic active center of lysozyme. The results of docking 
further showed that there was a hydrogen bond between 
Moby and lysozyme, so the interaction between Mobic and 
lysozyme was driven by hydrophobic interaction and 
hydrogen bond. 

Keywords — Mobic; lysozyme; spectroscopy; 

conformational; molecular docking. 

I. INTRODUCTION 

Non-steroidal anti-inflammatory dmgs can relieve pain and 
edema to play a role in the treatment of inflammation, 
however, patients taking non-steroidal anti-inflammatory 
dmgs often cause gastritis, gastric ulcer, kidney and liver 
damage and other adverse symptoms [1].Mobic is one of 
the non-steroidal anti-inflammatory dmgs [2] and the 
structural formula is shown in Figure 1. The toxic and side 

http://www.aipublications.com/iicmp/ 


effects of Mobic are much less than those of some other 
non-steroidal anti-inflammatory dmgs (such as dotaline, 
ibuprofen, etc.). Therefore, the adverse effect of Mobic on 
gastrointestinal function was much smaller [3] .This dmg is 
widely used in daily life and it is a very common dmg for 
the treatment of rheumatoid arthritis, in addition to the 
treatment of inflammatory diseases, Mobic also has a good 
effect on relieving physical pain in patients [4]. 



The relative molecular weight of lysozyme (LYSO) is 
about 14600 [5]. The main fluorescent groups in LYSO are 
Trp62 and Trpl08 amino acid residues [6].LYSO is widely 
found in humans and animals. Such as tears, saliva, blood 
and lymphoid tissue all contain LYSO, which is a very 
important antibacterial protein in organisms. In addition to 
antibacterial activity, LYSO also has the ability to treat 
inflammation, regulate the immune ability of the body, 
antihistamine and inhibit tumor activity [7]and many other 
biological functions. Because of its important physiological 
function, LYSO is often studied as an important protein 
model. 

In recent years, fluorescence spectroscopy has become 
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an important means to study the mechanism of 
ligand-protein system. So far, Sourav Das [8] has explored 
the relationship between naringin and LYSO by 
spectroscopy and molecular docking. Yan Juan [9]has 
studied the mechanism of interaction between vitexin and 
LYSO by fluorescence spectroscopy. ChenChen[10] has 
studied the binding of oxaliplatin to LYSO by fluorescence 
spectroscopy.At present, the methods to explore the 
mechanism of the reaction between drugs and proteins are 
mainly to investigate the fluorescence quenching of proteins 
at the excitation wavelengths of 280nmand 295nm with the 
change of drug concentration,which is to study the overall 
interaction between drugs and proteins [11]. And 
synchronous fluorescence has both good selectivity, high 
sensitivity and it can narrow band and reduce scattering 
interference! 12]. hi the study of the reaction mechanism 
between dmg and protein, synchronous fluorescence method 
is mainly used to investigate the conformational change of 
protein by the displacement of synchronous fluorescence 
peak [13].At present, there are few reports on the interaction 
between drugs and proteins by synchronous fluorescence 
spectrometry. In this paper, the interaction between Mobic 
and LYSO was studied by synchronous fluorescence 
spectroscopy, ultraviolet spectroscopy, circular dichroism 
spectroscopy and molecular docking simulation. The 
reaction mechanis m between drug molecules and amino acid 
residues can be revealed at a deeper level, which provides a 
certain reference for the study of the binding of small 
molecular drugs and proteins. 

II. EXPERIMENTAL 

2.1 Apparatus 

RF-5301PC fluorometer (Shimadzu, Japan); UV-3600 
UV-vis spectrophotometer (Shimadzu, Japan); 
MOS-450/SFM300 circular dichroism (Bio-Ix>gic, France); 
SYC-15 b super constant temperature water bath (Nanjing 
Sanli Electronic equipment Factory); SZ-93 automatic 
double Pure Water Distiller (Shanghai Yarong biochemical 
instrument Factory). 

2.2 Materials 

LYSO (purity grade inferior 99%, Sigma), reserve 
solution(5.0xl0' 6 mol/L); Mobic (CAS#,71125-28-7), 
reserve solution(4.0xl0" 4 mol/L).Tris-HCl buffer 
solutionwas used to keep the pH of the solution at 7.40, and 
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NaCl (0.10 mol/L) solution was used to maintain the ionic 
strength of the solution.The water used in the experiment 
was secondary quartz distilled water, and the above storage 
solutionwas kept away from light at 277 K.The fluorescence 
signal measured in the experiment was corrected by the 
''internal filter effect" Eq. (1)[14]: 

Fcor = Fobsxe (A - +Aem)n - 

( 1 ) 

Where F CO r and Fobs are the corrected and observed 
fluorescence signals, respectively, and A ex and A em are the 
absorbance values of Mobic-LYSO system at excitation and 
emission wavelengths, respectively. The fluorescence signal 
used in this article was corrected. 

2.3 Experiment procedure 

2.3.1 Synchronous Fluorescence measurements 

At 298K, 310K and 318K, 1.0 mL Tris-HCl buffer 
solution, 1.0 mL LYSO solution and different volume of 
Mobic solution were added to the 10.0 mL colorimetric tube 
at constant volume and constant temperature of 30 min. The 
prepared solution was placed in a 1.0 cm quartz colorimetric 
dish, and the difference A2 between the fixed emission 
wavelength and the excitation wavelength was set at 5 nm. 
The synchronous fluorescence spectra of the scanning 
system were 15 nm and 60 nm, respectively, and the 
synchronous fluorescence intensity was recorded. 

2.3.2 Circular dichroism measurements 

At 298K, 1.0 mL Tris-HCl buffer solution, 1.0 mL 
LYSO solution and different volume of Mobic solution were 
added to the 10.0 mL colorimetric tube at constant volume 
and constant temperature of 30 min.Circular dichroism 
measurements were performed with a 1.0 cm path length 
quartz cuvette. Each spectmmwas recorded at wavelengths 
between 190 and 300 nm and a scan speed of 1 nm/s. 

2.3.3 UV-Vis measurements 

At 298 K, 1.0 mL Tris-HCl buffer solution,1.0 mL 
LYSO solution and different volume Mobic solution was 
added to the 10.0 mL colorimetric tube at constant volume 
and constant temperature of 30 min. The absorbance of the 
system was determined by using the corresponding 
concentration of Mobic solution as the blank reference, and 
the UV absorption spectrum of the systemwas drawn. 

2.3.4 Molecular docking 

The crystal structure (PDB ID: 6LYZ) of LYSO comes 
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from the protein database (Protein Data Bank). The 
ChemDraw Pro 14.0 and ChemBio 3D Ultra 14.0 are used to 
draw the Mobic structure, and the energy minimization of 
the three-dimensional structure is carried out. AutoDock 


4.2.6 was used to study the molecular docking of Mobic and 
LYSO, and genetic algorithm was used to calculate the 
binding conformation of Mobic and LYSO[15]. 


III. RESULTS AND DISCUSSION 
3.1 Synchronous fluorescence studies of Mobic-LYSO system 




Fig.2: Synchronous fluorescence spectra of Mobic-LYSO system (T=298 K). (A) Al=60 nm; (B) AX=15nm 
Clyso=5.0x10' 7 mol/L, 1-7 C Mo bic=(0, 0.5, 1.0, 1.5, 2.0, 3.0, 4.0)xl0- 5 mol/L 


Synchronous fluorescence spectra can show the micro 
environmental information of protein fluorescence groups. 
When the A2 value between excitation wavelength and 
emission wavelength is stable at 15 nm or 60 nm 
synchronous fluorescence gives the characteristic 
information of Tyr residue or Trp residue, respectively 
[16] .As shown inFigure 2, the fluorescence of Tyr residue or 
Trp residue is quenched by the addition of Mobic, and the 
quenching degree is further strengthened with the increase of 
Table. 1: Quenching reactive parameters o 


drug concentration, which indicates that Mobic interacts 
with LYSO to form a stable complex.lt can be seen fromthe 
synchronous fluorescence spectra that the fluorescence 
peaks of Trp residues and Tyr residues in LYSO have shifted 
very weakly. This means that the binding of Mobic to LYSO 
changes the microenvironment of amino acid residues in 
small molecular proteins such as LYSO, which reduces its 
hydrophobicity and increases the extension of peptide chain 
[17]. 

bic-LYSO system at different temperatures 


A2 

(nm) 

77(K) 

A: t/ (L/mol-s) 

/Gv(Uinol) 

ri 

KaQJrmX) 

n 

r2 


298 

1.95xl0 12 

1.95X10 4 

0.9946 

2.13xl0 4 

1.04 

0.9907 

60 

310 

1.43xl0 12 

1.43X10 4 

0.9980 

1.58xl0 4 

1.05 

0.9952 


318 

1.23xl0 12 

1.23X10 4 

0.9936 

1.21xl0 4 

1.03 

0.9955 


298 

1.48xl0 12 

1.48X10 4 

0.9939 

1.51xl0 4 

1.06 

0.9951 

15 

310 

1.16xl0 12 

1.16X10 4 

0.9943 

1.13xl0 4 

1.09 

0.9956 


318 

0.91xl0 12 

0.91X10 4 

0.9926 

0.84xl0 4 

1.07 

0.9966 


n is the linear relative coefficient of Fo/F~[L \; n is the linear relative coefficient of \g[{Fo-F)l]~\g{[L\-n[Bt\{Fo-F)IFo} 


The Stem-\blmer [18] equation was shown below, through which the quenching constant Ks V and the quenching 
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rate constant k q can be calculated by using the fluorescence 
signal data obtained from the experiment: 

F 0 / F = 1 + k q r 0 m =1 + K SV IL\ ( 2 ) 

Where Fo and F represent the fluorescence signals in 
the absence and presence of quencher, respectively, to is the 
average lifetime of fluorescence without quencher, which is 
about 10" 8 s. Ksv is the Stem-Volmer quenching constant. k q 
is the bimolecular quenching constant, and [L] is the 
concentration of Mobic .According to the calculation of Eq. 
(2), the values of K sr and k q were Ms ted in Table 1. The results 
showed that the k q values at different temperatures are larger 
than the maximum diffusion collision quenching constant of 
2xl0 10 L/mol-s[19] for biomolecules by various quenching 
agents. At the same time, it could be seen from the data in 
Table 1 that with the increase of temperature, k q of 
Mobic-LYSO system and K sv 's are decreased. The results 
showed that the quenching mode of Mobic-LYSO system 
was static quenching. 

For static quenching, the Eq. (3)was generally used to 
calculate the binding constant K a and the number of binding 
sites ir. 


lg 


F n -F 


F, - F 


= n lg K a + n lg < [L\ - n — [5, ] \ (3) 


[Bt\ represents the concentration of LYSO, and the 
results were shown in Table l.From Table 1, n~l at the 
experimental temperature indicates that there is only one 
high affinity binding site[20] for Mobic to bind to LYSO, 
that is, Mobic formed a 1:1 complex with LYSO .The 
binding constant K a of Mobic to LYSO decreases with the 
increased of temperature, which further proved that the 
fluorescence quenching type of Mobic-LYSO system was 
static quenching.From the data in the table, it could be 
further seen that the K a value at A4=15 nm was obviously 
smaller than that at A2=60 nm, which indicated that the 
reaction between Mobic and Trp residues is stronger. 

3.2 Circular dichroism spectra studies of Mobic-LYSO 
system 

As one of the well-known biophysical techniques, 
circular dichroism spectrum (CD) is usually used to clarify 
the secondary stmcture of proteins, hi the CD spectrum, the 
far ultraviolet region (180 nm) mainly shows the «-helix, 
/(-fold and irregular crimping of the protein[21].The peaks 


at the left and right sides of the 208 nm and 222 nm are 
characteristic peaks of «-helix. The change of «-helix 
content can be obtained by using the following equation 
(4)[22]: 


K 


a-Helix 


(%) = ■ 


■N 


MRE 


-4000 


xlOO (4) 


33000-4000 

The Nmre is the molar ellipticity at 200 nm the 
number 4000 is the Nmre value of /(-fold and random curl at 
200 nm; and the number 33000 is the Nmre value of pure 
«-helix at 200 nm 


N 


CD 


MRE 


10 C P nl 


( 5 ) 


As shown in Eq. (5), Icd is the signal of CD at 208 nm, 
C P is the concentration of protein, n is the number of amino 
acid residues and / is the path length (cm). The circular 
dichroism spectroscopy of Mobic-LYSO system was shown 
inFigure 3. When the molar ratio of LYSO to Mobic were 
1:0, 1: 5 and 1: 10, the content of a-helix structure of LYSO 
molecules decreased from 22.76% to 10.79%. The intensity 
of the negative peaks decreases gradually, without obvious 
change in position and shape of the peaks. It could be 
concluded that the interaction between Mobic and LYSO 
made the a-helix stmcture of LYSO loose and changed the 
protein secondary stmcture, leading to quenching of the 
LYSO fluorescence. But a-helical stmcture was still 
dominant. 



Fig.3: The circular dichroism spectra of Mobic-LYSO 
system (T=298 K) 

Clyso= 2.0x10" 6 moI/L;CMobic= (0, 1.0, 2.0) xlO 5 moi/L 
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3.3 UV-vis absorption spectra studies of Mobic-LYSO 
system 



Fig.4: Absorption spectrum of Mobic-LYSO system (T=298 

K) 

C L Yso=2.0xl0 6 mol/L, 1-6: C M obic=(0, 0.2, 0.5, 1.0, 1.5, 
2.0)xl0 s moVL 

UV-vis absorption spectra can be used to explore the 
structural changes of proteins and to study the formation of 
protein-ligand complexes[23]. Figure4was an absorption 
spectrum of the Mobic-LYSO system From the figure, it 
can be seen that LYSO has two absorption peaks, and its 


strong absorption peak near 208 nm reflects the frame 
conformation of the protein. The weak absorption peak at 
about 280 nm is due to aromatic amino acids (Trp.Tyr and 
Phe)[24]. With the increases of Mobicconcentration, 
Figure4 showed that the intensity of the absorption peak at 
208nm decreased with the blue shift, and the absorption 
peak at 280nm also decreased slightly. This result indicated 
that the interaction between Mobic and LYSO led to the 
formation of new complexes, and the LYSO molecule 
tended to fold, the hydrophobicity of LYSO 
microenvironment was enhanced.The UV-Vis absorption 
spectrum of the system shifted, indicating that the 
fluorescence quenching of LYSO by Mobic, a new complex 
formed by the reaction, should belong to the static 
quenching process [25], which was consistent with the 
conclusion of synchronous fluorescence experiment. 

3.4 Type of interaction force of Mobic-LYSO system 

The thermodynamic parameters of Mobic-LYSO 
system we re calculated according to van't Hoff equation [26], 
and the calculated results were shown in Table 2. 

R\nK = AS-AH IT (6) 


AG = -RT In K-/AH-TAS 


(7) 


Table 2.The thermodynamic parameters of Mobic-LYSO at different temperatures 


System 

77(K) 

Ay(Umol) 

A///(kJ/mol) 

A5V(J/mol-K) 

AG/(kJ/mol) 


298 

2.13X10 4 


9.09 

-23.69 

AX= 

60 nm 

310 

1.58X10 4 

-21.98 

9.47 

-24.92 


318 

1.21X10 4 


9.03 

-24.85 

AX= 

298 

1.51X10 4 


3.90 

-23.84 

15 nm 

310 

1.13X10 4 

-22.68 

4.43 

-24.05 


318 

0.84X10 4 


3.81 

-23.89 


Where R is a gas constant (AH and AS with a value of 
about 8.314). Mobic and LYSO could be calculated by a 
linear relationship between the natural logarithm fin K„) of 
the binding constant and the reciprocal (1/T) of the 
temperature. The results were shown in Table 2.It could be 
seen from Table 2 that AG<0 indicated that the binding 
reaction between Mobic and LYSOwas spontaneous, and 
AH<0 indicated that the formation of Mobic-LYSO complex 
was exothermic .The arrangement of water molecules created 
a more random configuration around drugs and proteins in a 


more orderly manner. Therefore, A.S’>0 was usually used as 
evidence of hydrophobic interaction between dmg 
molecules and protein molecules[27]. Based on this, the 
hydrophobic interaction between Mobic and LYSO could be 
judged. 

Ross and Subramanian [28]believed that when AH~0, 
AS>0, there was electrostatic attraction between dmg 
molecules and biomolecules, but now some reports thought 
that when AH<0, AS>0 could directly judge that the main 
type of force between the binding systems was electrostatic 
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interaction [29]. In order to further verify whether the main 
force between Mobic and LYSOwas electrostatic force, the 
effect of ionic strength on Mobic-LYSO interaction was 
discussed in this paper. 

The concentrations of Mobic and LYSO were fixed, 
and different concentrations of NaCl were added, and F/Fo 
was plotted against CNaa and the results were shown in 
Figure5. The experimental results showed that when the 
concentration of NaCl increases, the ratio of F/Fo did not 
change significantly. The results showed that the binding of 
Mobic to LYSOwas not affected by ionic strength, that was, 
the electrostatic interaction of Mobic-LYSO system was not 
obvious. If the electrostatic interaction played a leading role 
in the binding of protein to ligands, as the increasing of the 
salt concentration, the interaction intensity between protein 
and ligand decreases gradually [30]. This also indicated that 
when AH<0, AS>0 could not directly judge the main force 
between Mobic-LYSO system was electrostatic force. 



Fig.5: Synchronous fluorescence intensity ofMobic-LYSO 
system as a function of NaCl concentration (T=298 K) 


Clyso=5.0x10' 7 mol/L, CMobic= 1.0x1 O' 5 mol/L, CNaci= (0, 
0.1, 0.2, 0.5, 1.0, 1.5, 2.0, 2.5, 3.0) xlO 1 mol/L 

3.5 Molecular docking 

Molecular docking plays an important role in 
exploring the interaction between ligands and receptois. In 
order to further determine the binding position of 
Mobic-LYSO system and the effect of Mobic binding to 
LYSO on ligands and receptois. In this paper, the binding 
model of Mobic and LYSOwas established by molecular 
docking method. By this method, the type of force and the 
lowest binding energy of Mobic and LYSO binding system 
could be obtained.Figure6(A) showed the optimal binding 
position of Mobic to LYSO, where the Trp62 residue and 
the Trp63 residue formed two hydrogen bonds with Mobic, 
respectively, with bond lengths of 2.105 A and 1.809 A, 
respectively. This result showed that hydrogen bond played 
an important role in the binding of Mobic to LYSO. 
Figure6(B) showed a plurality of hydrophobic amino acid 
residues such asAla95,Alal07, Val99 and Vall09around 
Mobic, further indicating that there was hydrophobic force 
in the binding process between Mobic and LYSO.The 
amino acid residues such as Tyr53, Trp62andTrp63 were 
relatively close to the binding position of Mobic and LYSO, 
which led to the binding could effectively quench the 
endogenous fluorescence of LYSO, which was consistent 
with the conclusion of fluorescence quenching experiment. 
Glu and Aspare the key residues of the catalytic active 
center of LYSO[31]. The results of molecular docking also 
showed that the binding of Mobic and LYSO could change 
the microenvironment of the catalytic active center of 
LYSO.In other words, the binding of the system might 
affect the catalytic activity of LYSO. 



Fig.6: Molecular docking of the Interaction between Mobic and LYSO. 

(A) Mobic located within the hydrophobic pocket in LYSO(B) Detailed illustration of the amino acid residues lining the 

binding site in theMobic and LYSO cavity 
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The binding energy obtained from molecular docking 
for Mobic and LYSO interaction was -26.8lkJ/mol. At 
310K, it is close to the synchronous fluorescence 
experiment that AG (Trp)=24.92 kJ/mol, AG (Tyr) =24.05 
kJ/mol.This difference may be due to exclusion of the 
solvent in docking simulations or rigidity of the receptor 
other than Trp and Tyr residues [32]. The energy data 
obtained by docking the molecules were list in Table 3. 
From Table 3, it could be also seen that the electrostatic 
energy was very much lower than the sum of van der Waals 


energy, hydrogen bonding energy and desolvation free 
energy in the binding process of Mobic with LYSO, 
indicating that the main interaction mode between Mobic 
and LYSOwas not electrostatic binding mode.Combined 
with the data of fluorescence experiments and the results of 
theoretical modeling, it could be seen that hydrophobic 
interaction and hydrogen bond were the main forces driving 
the combination of Mobic molecules with LYSO molecules, 
which led to the static quenching of LYSO. 


Table 3. Docking energy ofMobic-LYSO system (unit: kJ/mol) 


Protein PDB ID 

A Go 

A Ei 

A E 2 

A E 3 

6LYZ 

-26.81 

-30.53 

-30.32 

-0.21 


A Go is the binding energy in the binding process. 

A Ei denotes intennolecular interaction energy, which is a sum of van der Waals energy, hydrogen bonding energy, 


desolvation free energy and electrostatic energy. 

AE:is the sum of van der Waals energy, hydrogen bonding energy and desolvation free energy. 
A£jis the electrostatic energy. 


IV. CONCLUSION 

Under simulated physiological conditions, a variety of 
spectroscopy and molecular docking techniques were used 
to explore the interaction between Mobic and LYSO, which 
is ofpositive significance forthe study of pharmacokinetics, 
pharmacodynamics and toxicology of Mobic.At present, the 
use of spectroscopy to study the interaction between proteins 
and ligands is more and more popular and favored by 
researchers. Compared with other research methods such as 
chromatography, this research method has the advantages of 
low experimental cost and short measurement period. 
Therefore, this study provides an important reference forthe 
efficient and rapid study of the interaction between proteins 
and ligands. 
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